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1 Trapped Bose atoms 



Trapped gases of Bose atoms constitute an important class of nonuniform systems. Ther- 
modynamics and dynamics of trapped atoms have been studied in many details. For review, 
see the book [1] and review articles [2-5]. One of the topics that have met a great deal 
of controversy is the problem of particle fluctuations in the systems with Bose-Einstein 
condensate. There have appeared a number of publications claiming the occurrence of 
thermodynamically anomalous condensate fluctuations everywhere below the condensation 
temperature, for both uniform as well as trapped systems. Thermodynamically anoma- 
lous fluctuations correspond to the particle dispersion proportional to iV 4 / 3 , instead of 
N for the normal fluctuations. As has been explained in the review papers [3,5], the 
occurrence of these thermodynamically anomalous fluctuations implies thermodynamic in- 
stability. Hence, the systems with such anomalous fluctuations simply cannot exist. For 
the case of uniform systems, it has been shown [6-8] that the thermodynamically anoma- 
lous fluctuations are due to incorrect calculations, while the correct calculational procedure 
yields thermodynamically normal particle fluctuations. 

In the present paper, this result is generalized to the case of arbitrary nonuniform sys- 
tems, including the case of trapped atoms. It will be shown that for any Bose-condensed 
system, whether uniform or nonuniform, particle fluctuations are always thermodynami- 
cally normal. The notion of thermodynamically normal fluctuations is related to that of 
thermodynamic limit. For a uniform system of N atoms in volume V, the thermodynamic 
limit is commonly defined as 

N 

N oo , V — > oo , — — > const . (1) 

For an arbitrary nonuniform system, the thermodynamic limit can be defined [9] as the 
limit 

(A) 

N -> oo , (A) -> oo , -> const , (2) 

valid for the statistical average < A > of any extensive observable A. 
Atoms are often trapped by means of a power-law confining potential 



2 

Q=l 



(3) 



where d is space dimensionality, n a > 0, and the potential parameters are related by the 
equations 

mli ^mu a 

It is convenient to introduce the effective frequency and the effective localization length, 
respectively, 

/ d \ l ' d 

z °=(n z °) ' ® 

which are connected by the relations 

1 , 1 
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For this type of the confining potentials, the thermodynamic limit (2) reduces [9] to 
the limit 

iV — > oo , ujq — > , Nuq — > const , (6) 

or, equivalently, to 



N 
P 



N — > oo , Zq^oo, TU - ^ const , (7) 



where the notation of the confining strength 
is introduced. 

The passage to the uniform system corresponds to n a — > oo, when s — > d/2 and 
/o — > £/2, with L being the system length, such that L d = V. 

The comparison of Eqs. (1) and (7) tells us that, for a confined system, the effective 
volume can be defined [1,9] as V = const -II s . This definition seems to be not unique, since 
the proportionality constant here is not yet prescribed. However, we can remember that 
real trapped systems are always finite, being bounded by the trap volume V . The power- 
law confining potential (3) is just a model for the real trapping potential. For sufficiently 
large traps, this is a good model, which does not contradict the fact that the real trap has 
a finite volume V. Therefore, the relation V = const ■ /q s can be treated as defining the 
proportionality coefficient. 

If H is the system grand Hamiltonian, then the grand thermodynamic potential is 

n = -PV = -T In Tr e~ pH , (9) 
where T is temperature and j3 = 1/T. This defines the system pressure 

P = 1 p In Tr e-? H , (10) 

in which p = N/V is the average density of atoms. Thus, remembering that any trap has 
a finite volume V allows us to use, for trapped atoms, the same thermodynamic relations 
as for uniform systems. 

Particle fluctuations are characterized by the dispersion 

A 2 (N) = (N 2 ) - (N) 2 (11) 

of the number-of-particle operator^. This dispersion is straightforwardly connected with 
the isothermal compressibility 

1 fdPY 1 A 2 (N) , x 

hydrodynamic sound velocity 

1 (dP\ 1 NT 

s T = — — = = ^ , 13 

m\dp J mpK T mA 2 (N) 
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where m is atomic mass, and with the central structure factor 



T A 2 (N) 

s (°) = T ^=^4 = ^ i ' < 14 > 

Equations (12), (13), and (14) are exact, being valid for any nonuniform system. In order 
that these measurable quantities would have sense, it is necessary and sufficient that the 
particle fluctuations be thermodynamically normal, such that 

^ ' -> const (N -> oo) . (15) 

If the particle fluctuations would be thermodynamically anomalous, being proportional to 
iV 4 / 3 , as is claimed by some authors, the compressibility (12) would be divergent, sound 
velocity (13) would be zero, while the structure factor (14) would be infinite. Such a 
thermodynamically anomalous behavior would mean that the considered system is ther- 
modynamically unstable. 



2 Nonuniform Bose systems 

An arbitrary nonuniform Bose system is characterized by the energy Hamiltonian 

H = J ft(r) ^_ + uj V>(r) dr + 

+ ^ J ^ t (r)^ t (r / )$(r - r')V>(r')V>(r) drdr' , (16) 

in which ■?/> is the Bose field operator, U — U(r) is an external potential, such as the 
trapping potential, and $(r) is an integrable symmetric interaction potential. Here and in 
what follows, the units are used where h — 1, ks — 1. 

The appearance of the Bose-Einstein condensate is the necessary and sufficient condition 
for the global gauge symmetry breaking [10]. The latter is conveniently realized by means 
of the Bogolubov shift [11, 12] for the field operator 

^(r)=77(r)+^(r). (17) 

Here 77 (r) is the condensate wave function and ipiij) is the Bose field operator of uncon- 
densed atoms. To exclude the double counting, these variables are to be orthogonal to 
each other: 

V(r)V>i(r) dr = . (18) 



I 



The condensate wave function is normalized to the number of condensed atoms 

N = J p (r)dr, po(r) = |r/(r)| 2 . (19) 
And the number of uncondensed atoms is given by the average 

N, = (A\) , N ± = J 4(r)^(r) dr . (20) 
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The total number of atoms in the system is 

N = N + JVi = (N) , N = N + jVi . (21) 

The condensate wave function plays the role of the order parameter characterizing the 
gauge symmetry breaking: 

77(r) = $(r)) , (^i(r)) = . (22) 

This means that the Hamiltonian should not contain the terms linear in the field operators 
of uncondensed atoms [13]. The latter condition can be realized by complimenting the 
Hamiltonian with a counterterm 



A 

for which 



J [A(r)V>!(r)+A*(r)V>i(r)] dr , (23) 



(A) = , (24) 

and the Lagrange multipliers A(r) are chosen so that to kill the terms linear in ipi(r). 

Taking into account the statistical constraints (19), (20), and (24) defines the grand 
Hamiltonian 

H[r], fc] = H- fi N - ^N, - A , (25) 

which is a functional of the field variables rj and ipi- The quantities /i and /ii are the 
Lagrange multipliers guaranteeing the validity of the normalization conditions (19) and 
(20). With the Bogolubov shift (17), Hamiltonian (25) is the sum 

4 

71=0 

of the terms labelled according to the entering powers of the operators ipi. The zero-order 
term 

if <°) = J V *(r) (- ^ + U - fj^j 77(r) dr + 

+ l - J $(r - r') \r](r)\ 2 ^(r')) 2 drdr' (27) 

does not contain the field operators of uncondensed atoms. The first-order term ff W = 0, 
being eliminated by the linear killer (23). The second-order term is 

H (2) = J (- ^ + U - ^ Vi(r) dr + 

)] (r)|V 1 t (rVi(r')+^(r)^r , )^ 1 t (r')^i(r)+ 
+ l - rf (r)?f (r>i(r>i(r) + I 7 ? (r)r ? (r / )^ 1 t (r')^ 1 t (r)l drdr' . (28) 



+ / $(r-r' 
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Respectively, the third-order term is 

= y$(r-r') ?7*(r)V'I(r / )V'i(r / )V'i(r)+V'i t (r)V'i t (r / )V'i(r / )'7( r ) d *dr' , (29) 
and for the fourth-order term, one has 

#( 4 ) = 1 y v4(r)V4(r')$(r - rOVitrOV'iM drdr' . (30) 

The evolution equations for the variables rj and ipi can be represented by the variational 
forms 

and 

z — ^i(r,t) = (32) 
#j(r,t) 

Let us introduce the normal density matrix 

Pl (r,r') = {4(r')Mr)) (33) 

and the so-called anomalous density matrix 

ax(ry) = (^(rVxO:)) . (34) 

Their diagonal elements define the density of uncondensed atoms and the anomalous av- 
erage, respectively: 

pi(r) = (Vi f (r)Vi(r)) , <7i(r) = (^i(r)^(r)) . (35) 

The local density of atoms is the sum 

p(r)=p (r)+ Pl (r). (36) 

We shall also need the triple correlator 

e(r,rO = (^(r')^(r')^(r)). (37) 

With these notations, Eq. (31) results in the evolution equation for the condensate 
wave function 

d / V 2 \ 
i — r](r) = - \-U-no] r](r) + 



dt w V 2m 

+ J $(r - r') [p(r')77(r) + Pl (r, r'Mr') + ^(r, r')^(r') + £(r, r')] dv' , (38) 
while Eq. (32) yields the equation of motion for the operator of uncondensed atoms 



l5Wr)=(~ + "-MW') + 



+ J $(r-r') ^(r')| 2 ^i(r) + 77*^)77 (r)^i(r') + ^(^(r^iV) + X(r,r')J dr' , (39) 
in which 

X(r,r') = 4(r')Mr'Mr) + 4(r')v(r')M^) + V* + ^(r>i(r>i(r) ■ 
For an equilibrium system, one has 

d 

— r)(r) = (equilibrium) . 
Then Eq. (38) reduces to the eigenvalue problem 

"»" (r) = (- S + u ) " (r) + 

+ j *(r - r') [p(r>(r) + Pl (r, r')i)(r') + ^(r, r')i(*(r') + £(r, r')] dr' , (40) 

defining the condensate wave function and the Lagrange multiplier fiQ. 

Note that the Lagrange multipliers /io and do not need to coincide with each other. 
Their relation to the system chemical potential /i is given [13-15] by the equality 

fi = fi n + /Jini , 

in which the atomic fractions of condensed and uncondensed atoms, respectively, are 

The formalism of this section provides the basis for the self-consistent theory of arbitrary 
nonuniform Bose-condensed systems [8,13-15]. 



3 Hartree-Fock-Bogolubov approximation 

To proceed further, let us resort to the Hartree-Fock-Bogolubov (HFB) approximation, 
following the way of Refs. [13-15]. Then Hamiltonian (26) in the HFB approximation 
becomes 

H HFB = E HFB + J ^J(r) (- ^ + U - /i^j ^i(r) dr + 



+ J $(r - r') 



p(r')W(r)^i(r) +p(r',r)^(r')^(r) + - a(r, r')i/j\(r')4(r)+ 



drdr' , 



(41) 



where 



E HBF = H^- l - J $(r-r') [ Pl (r) Pl (r') + | Pl (r, r')| 2 + (^(r, r')| 2 ] drdr' , (42) 



7 



and the notation is introduced for the total normal density matrix 

p(r,r) = r](r)r]*(r')+ Pl (r,r') > 
and for the total anomalous density matrix 

cr(r, r') = rj(r)r](r') + <7i(r, r') . 
In the HFB approximation, the condensate-function equation (38) takes the form 



(43) 



(44) 



1 at ^ = 



2m 



+ U-fi ) r](r) + 



+ J <t>(r - r') [p(r')rj(r) + Pl (r, r>) V (r') + ^(r, r') V *(r')} dr' . 
And Eq. (39) for the field operator of uncondensed atoms reduces to 

V 2 



(45) 



2m 



+ U-pL 1 ) ^(r) + 



/ 



+ /$(r-r') p(r')Vi(r)+p(r, r'^r') + a(r,r')4(r') 



dr' . 



(46) 



For an equilibrium system, the eigenvalue problem (40) becomes 



+ 



V 2 \ 
- - + Uj ,(,) + 

/*(r - r') |p(r')r,(r) + Pl (r, 1>|Y) + a,(r, r')r,'(r')] dr' . (47) 

The HFB Hamiltonian (41) can be diagonalized by means of the general Bogolubov 
canonical transformations [16] that, in the used notation, read as 

^i( r ) = ^[W r ) + 6 k(r)] , 

k 

<(r)Vi(r) - T£(r)irf(r)] dr' , (48) 
where A; is a multi-index. Since ipi is the Bose operator, one has 

£ MrK(r') - ^(rH(r')] = *(r - r') , 

^K(rK(r')-^(rK(r')] = 0. (49) 

fc 

And requiring that be also a Bose operator gives 
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/ 



J [u k (r)v p (r) - v k (r)u p (r)] dr = . 



(50) 



The condition that Hamiltonian (41) be diagonalized by the Bogolubov transformation 
(48) is equivalent to the Bogolubov equations 

J Kr,r'K(r') + A(r, r> fe (r')] dr' = e k u k (r) , 
j [u*(r, r')v k (r') + A*(r, r')u k (r')} dr' = -e k v k (r) , (51) 



in which 



u(r,r') 



^ + U(r) - ^ + J $(r - r')p(r') dr' 



5(r - r') + 



+ $(r - r')p(r, r') (52) 

and 

A(r,r') = $(r-r>(r,r') . (53) 

The Bogolubov equations (51) is the eigenproblem defining the Bogolubov functions u k (r) 
and v k (r) and the spectrum of collective excitations e k . 

As a result of the diagonalization, Hamiltonian (41) reduces to the Bogolubov form 

H B = E B + J2e k b{b k , (54) 
k 

where 

E B = E HFB ~J2 £k J M r )! 2 dr ■ ( 55 ) 

k 

With the diagonal Hamiltonian (54), it is straightforward to calculate the required averages. 
Thus, for the distribution of excitations, one gets 

7 r fc ^(^) = (e^-l)- 1 . (56) 
The normal density matrix (33) is 

Pl (r, r') = WkM'W) + (1 + *kK(r)v k (r')} . (57) 
k 

And the anomalous density matrix (34) becomes 

<7!(r, r') = [KkU k (r)v* k (r') + (1 + 7r k )v* k (r)u k (r')} . (58) 
k 

For the diagonal elements of Eqs. (57) and (58), we obtain 

pi(r) = [^\u k (r)\ 2 + (1 + n k )\v k (r)\ 2 ] (59) 

k 

and, respectively, 

( 7 1 (r) = ^(l + 27r fe K(rK(r). (60) 
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4 Local-density approximation 

When the spatial variation of the external nonuniform potential is sufficiently slow and the 
trap is sufficiently large, one can employ the local- density approximation [1,2,17], which 
is also called the semi-classical approximation and is widely used for describing trapped 
atoms [9,18,19]. 

In this approximation, one looks for the Bogolubov functions represented as the prod- 
ucts 

u fc (r) = u(k, r)(p k (r) , v k (r) = v(k, r)y? fc (r) (61) 
factorized with the plane waves 

Mr) = ^ ■ (62) 



and 



The amplitudes u(k, r) and v(k, r) are assumed to be slowly varying in space, as compared 
to the spatial variation of the plane wave, 

|V^(k,r)| |Vy fc (r)| ^ 
Kk,r)| ^ |^(r)| • 

In the Bogolubov equations (51), one makes the replacements 

J u(r,r')u k (r' ) dr' = u(k, r)u k (r) , 

J u(r, r')v k (r') dr' = cu(k, r)v k (r) , (64) 
J A(r,r'K(r') dr' = A(r)u k (r) , 

J A(r,r'K(r')rfr' = A(rK(r), (65) 

u(k, r) = + U(r) + 2$ oP (r) - ^(r) (66) 
is used, instead of Eq. (52), and the quantity 

A(r) = bo(r) + <7i(r)]$ (67) 
is used instead of Eq. (53), with the notation 

> = J $(r) dr . (68) 



where the form 



In what follows, we assume that $o > 0. Then the Bogolubov equations (51) reduce to the 
eigenproblem 

[u(k, r) - e(k, r)] u(k, r) + A(r)u(k, r) = , 
A*(r)u(k, r) + [co*(k, r) + e(k, r)] v(k, r) = . (69) 
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The amplitudes w(k, r) and i>(k, r) can be taken to be real and, in view of Eqs. (49), 
satisfying the relation 

w 2 (k,r) -v 2 (k,r) = 1 . 
Solving eigenproblem (69) yields the Bogolubov spectrum of collective excitations 

e(k, r) = v/^ 2 (k,r)-A 2 (k,r) . (70) 

And for the amplitudes, we find 

2/i \ _ ^( k > r ) + ^( k > r ) 2/i , _ ^(Kr) - e(k,r) 



u 2 (k,r) = v V/i > « ( k > r ) 



2e(k,r) ' v ' ' 2e(k,r) 

„.(k,r W k,r) = -^L. (71) 

The necessary condition for the condensate existence [5], equivalent to the Hugenholtz- 
Pines theorem [20], requires that 

lim e(k, r) = , e(k,r)>0. (72) 

k— >0 

This gives the Lagrange multiplier 

Mr) = C/(r) + [p (r) + 2 Pl (r) - ^(r)]^ . (73) 
Introducing the local sound velocity c(r) by the equation 

mc 2 ^ [,oo(r) + <7i(r)]$ (74) 

transforms Eq. (66) to 

k 2 

w(k,r) = - + mC 2 (r), (75) 

while Eq. (67) becomes 

A(r) = mc 2 (r) . (76) 
For the Bogolubov spectrum (70), we obtain the familiar expression 



5(k,r) = yc 2 (r)P+^ 2 , (77) 

but with the local sound velocity. 

Instead of distribution (56) for excitations, we have 

7r(k,r) = {exp[/3e(k,r)]-l}- 1 , (78) 

with the symmetry properties 

7r(-k, r) = 7r(k, r) , e(-k, r) = e(k, r) . (79) 

Using Eqs. (71), we find the normal density matrix (57) as 

p 1 (r,r') = ^n(k,rK(rM,(r / ), (80) 

k 

11 



with the distribution of atoms 



. uCk, r) 
n(k, r) = — r coth 



g(k, r) 
2T 



2e(k,r) 

For the anomalous density matrix (58), we get 

<7i(r,r') = ^V(k,rV fc (r)^(r') , 



(81) 



(82) 



where 



a(k,r) = - 



mc 2 (r) 



coth 



g(k, r) 
2T 



(83) 



2e(k,r) 

The diagonal elements of Eqs. (80) and (82) give the density of uncondensed atoms and 
the anomalous average, respectively: 



(84) 



The grand thermodynamic potential (9) becomes 

n = E B + T J In [1 - exp{-/?5(k, r)}] 



rfk 



(85) 



Employing the above formulas, one can calculate all thermodynamic characteristics. 



5 Particle fluctuations and stability 

Particle fluctuations are defined by the number-of-particle operator dispersion (11). The 
latter is proportional to the isothermal compressibility (12) which is to be thermodynami- 
cally normal for the stability of the system. 

In the grand canonical ensemble, used here, the compressibility is 



Kt 



A 2 (N) 
pTN 



P = y I P(r) dr 



(86) 



In the canonical ensemble, the number of particles is fixed. However, it would be incorrect 
to conclude that the compressibility then is zero. Expression (86) is not defined for the 
canonical ensemble. Instead, one should use the formula 



K T = 77 



I ( d2F 

V [W 2 



TN 



where F is the canonical free energy. 

In the same way, it would be principally wrong to state that, fixing the number of atoms 
iV and setting the particle dispersion (11) to zero, would result in the equality A 2 (JV ) = 
A 2 (A r 1 ) that would define the condensate fluctuations A 2 (iVo) in the canonical ensemble by 
calculating the dispersion A^iVx) of uncondensed atoms. This way of reasoning is wrong 
because the particle dispersion (11) is not defined for the canonical ensemble. In addition, 
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calculating A 2 (iV"i) in the Fock space has nothing to do with the canonical ensemble, as 
far as in the Fock space, the number of particles is not fixed. 

In the Fock space, the correct conclusion [10] for the particle fluctuations, after using 
the Bogolubov shift (17), is that the condensate fluctuations are zero, 

A 2 (iV ) = , (87) 

and that the particle fluctuations are completely due to those of uncondensed atoms, 

A 2 (N) = A 2 ^) . (88) 

Expressions (87) and (88), for large particle numbers N are asymptotically exact [10]. 

It is also important that the fluctuations of the total number of particles are thermody- 
namically normal if and only if the fluctuations of both the condensate fraction as well as 
of uncondensed atoms are thermodynamically normal. And, vice versa, the fluctuations of 
the total number of particles are thermodynamically anomalous if and only if a least one 
of the fractions produces thermodynamically anomalous fluctuations. The corresponding 
theorem has been rigorously proved in Refs. [6-8]. 

When the HFB approximation is involved, Eq. (88) is not convenient to use, since the 
HFB approximation is an effective second-order approximation with regard to the powers 
of the field operators of uncondensed atoms ipi. In that second-order approximation, the 
quantity Nf is not defined, being of the fourth order with respect to ipi. But we can employ 
another way of calculating the particle dispersion (11). We may notice that the latter can 
be expressed through the density-density correlation function 

£>(r,r') = (ft(r)iP(r)ft(r')4>(r')) (89) 

as 

A 2 (iV) = J [D(r, r') - p(r)p(r')] drdr' . (90) 

The density-density correlation function (89) is related to the diagonal correlation function 

C(r,r') = (^(r)^( r ^( r ^( r )) (91) 

by the equality 

D(r, r') = p (t)6(t - r') + C(r, r') . (92) 

In turn, the diagonal correlation function (91) is connected with the pair correlation func- 
tion 

_ (jt(r)£(r3j(r0j(r» 
)_ PMp(r') (93) 

through the relation 

C(r,r')=p(r)p(r')g(r,r') . (94) 
Therefore the particle dispersion (90) can be represented as 

A 2 (iV) = N + J p(r)p(r') [g(r, r') - 1] drdr' = 
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N + J [C(r, r') - p(r)p(r')] drdr' 



(95) 



This formula is valid for arbitrary nonuniform systems. 

Accomplishing the Bogolubov shift (17) in Eq. (91) yields 



C(r,r')=p(r)p (r')+p (r)pi(rO + 

+2Re [77*(r)77*(r>i(r, r') + ?7*(r)?7*(r>i( r , r) + ?f (r)£(r, r') + ?f (r)£(r', r)] + 

+ C 1 (r,r'), (96) 

where 

d(r,r') = (^(0^(0^(0^1 W) ■ (97) 

In the HFB approximation, the triple correlator (37) is zero, while the correlation function 
^97) becomes 

diry) = Pl (r) Pl (r') + | Pl (r,r')| 2 + ki(r,r')| 2 . (98) 
Then the correlation function (96) reduces to 

C(r, r') = p{v)p{v') + 2Re [ V * (r) V (r') Pl (r , r') + ^(r)^ (rV^r, r')] + 

+ | Pl (r,r')| 2 + Mr,r')| 2 . (99) 

As is stressed above, the HFB approximation is of second order with respect to the oper- 
ators of uncondensed atoms. The terms, containing higher orders are not defined in this 
approximation and have to be omitted. This concerns the last two terms in Eq. (99). At 
the same time, in the frame of this approximation for an equilibrium system, the function 
r/(r) in the second and third terms of Eq. (99) can be replaced by the real value a/ p(r). 
As a result, for the particle dispersion (95), we obtain 

A 2 (A) = N + 2 J Vp(r)py) [Pi(r, r') + ^(r, r')] drdr' . (100) 

In the spirit of the local-density approximation, the latter expression can be rewritten as 

A 2 (A) = N + 2 J p(r) [ Pl (r, r') + ^(r, r')] drdr' . (101) 

Invoking Eqs. (80) and (82) gives 

/ pi(r, r') dr' = lim n(k, r) , / <7i(r, r) dr' = lim a(k, r) . (102) 

/ k — ^0 / k — ^0 

Hence dispersion (101) reads as 

A 2 (A) = A + 2 / p(r) lim [n(k, r) + a(k, r)} dr . (103) 



Using Eqs. (81) and (83), we get 



lim [n(k, r) + <j(k, r)] = - 

k— >0 Z 



1 - i 



mc 2 (r 
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It is worth emphasizing that both Eqs. (81), as well as (83), diverge in the long- wave limit 
k — > 0. But their divergences, being of opposite signs, cancel each other. This means that 
taking into account the anomalous average (83) is of principal importance. Without it, 
the dispersion (103) would diverge, and the compressibility (86) would be infinite, which 
implies the system instability. 

In that way, the particle fluctuations are described by the dispersion 



These formulas show that particle fluctuations are thermodynamically normal and com- 
pressibility (105) is finite, even in the thermodynamic limit. At the critical temperature 
T c , where c(r) tends to zero, the compressibility diverges as T — > T c . This divergence 
is typical for the point of a second-order phase transition, where the system is unstable. 
However everywhere below T c , the compressibility is finite and fluctuations are thermody- 
namically normal. This conclusion does not depend on the statistical ensemble involved. 
Both canonical and grand canonical ensembles for Bose systems produce the results co- 
inciding in the thermodynamic limit [21], provided all calculations are done correctly, by 
using the representative statistical ensembles [8,13-15,22]. 

The compressibility (105) is a directly observable quantity, being related to the structure 
factor (14). The latter can be measured by studying light scattering from ultracold atomic 
gases [23,24]. No one scattering experiment with Bose-condensed systems has ever showed 
a thermodynamically anomalous structural factor. 

Thermodynamically anomalous fluctuations for a stable statistical system can arise 
solely because of incorrect calculations. For example, if in the HFB or Bogolubov ap- 
proximation, one includes in calculations the fourth- or higher-order terms, which are not 
defined in the second-order approximation, then one can get any kind of thermodynamic 
anomalies. However, such anomalies are physically meaningless, being mathematically 
wrong. For correctly calculating the fourth-order terms, one has to use a fourth-order 
Hamiltonian. 

Nonperturbative thermodynamics of an interacting Bose gas, for all temperatures below 
T c , has recently been studied by Floerchinger and Wetterich [25] by using renormalizat ion- 
group techniques, which effectively take into account all orders of field operators. Their 
results confirm that the compressibility is finite everywhere below T c , hence, all particle 
fluctuations are thermodynamically normal. 
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(104) 



Therefore, the compressibility (86) is 




(105) 
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